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Abstract 

The structural and electronic properties of small scandium oxide clusters ScnOm (n=l-3, m=l- 
2n) are systematically studied within the screened hybrid density functional theory. It is found 
that the ground states of these scandium oxide clusters can be obtained by the sequential oxidation 
of small "core" scandium clusters. The fragmentation analysis demonstrates that the ScO, SC2O2, 
SC2O3, SC3O3, and SC3O4 clusters are especially stable. Strong hybridizations between 0-2p and 
Sc-Sd orbitals are found to be the most significant character around the Fermi level. In comparison 
with standard density functional theory calculations, we find that the screened hybrid density 
functional theory can correct the wrong symmetries, and yield more precise description for the 
localized 2>d electronic states of scandium. 

PACS numbers: 73.22.-f, 36.40. Cg, 36.40. Qv, 71.15.-m. 
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I. INTRODUCTION 



Transition metal oxide clusters have attracted enormous attention because their struc- 
tural, electronic, and magnetic properties are often quite different from those in their bulk 
phasei"— . For example, small clusters like (ZnO)„, (V205)„, (Cr03)„ form planar structures 
with very small sizes, and some small clusters present novel magnetic properties'*^"—. For 
scandium (Sc), continuous interests are shown for several reasons. Firstly, Sc is the first 
transition metal in the periodic table, and thus is always taken as a prototype to study the 
complex phenomena associated with d shell electrons^^. Secondly, Sc oxide nanostructures 
can be used as catalysts in some reactions like the selective reduction of nitric oxides with 
methanei^i^. Thirdly, some Sc oxide clusters have been recognized in the spectra of M-type 
stars^^ii^. Finally, some Sc oxide clusters can be steadily encapsulated into closed carbon 
cages of fullerenes to form a novel nanostructure^""— . 

The small Scandium oxide clusters can be prepared by laser ablation of scandium metal 
in the presence of oxygen-saturated atmosphere^"—. With the presence of oxygen or ni- 
trogen oxide gases, scandium oxide clusters ranging from ScO to ScsOg have already been 
generated and detected^^. Many experimental studies have already been applied to study 
the structural, energetic, vibrational, electronic, and magnetic properties of scandium oxide 
clusters. For example, the photoionization spectra of Sc„0 were measured and strongly size- 
dependent ionization potentials were observed22i2S. For the ScO molecule, the electron-spin 
resonance and optical spectroscopy in neon and argon matrices revealed that its ground state 
is a doublet^i^, and the molecular bond length and vibrational frequency have been mea- 
sured to be 1.668 A and 965 cm^^ respectively^^i^. Accompanying with these experimental 
results, lots of theoretical calculations, especially those based on density functional theory 
(DFT), were also carried out to explore the ground-state electronic structures of scandium 
oxide clusters. However, most of those studies on Sc oxide clusters were mostly on the 
standard DFT level, within the local density approximation (LDA) or generalized gradient 
approximation (GGA). In this paper, we systematically study the atomic and electronic 
structures of Sc oxide clusters using the screened hybrid density functional theory, which 
has proven to be able to significantly improve the description of finite, molecular systems, 
such as atomization energies, bond lengths and vibrational frequencies using standard DFT 
theories^^. 
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II. COMPUTATIONAL METHOD 



For extended systems, the ground-state electronic properties can be obtained by solving 
the Kohn-Sham equation within DFT— i^, utilizing the standard approximations for the 
exchange-correlation (xc) energy functional E^c, i-e., LDA and GGA. One problem of the 
two most commonly applied functionals is that they rely on the xc energy per particle of 
the uniform electron gas, and thus are expected to be useful only for systems with slowly 
varying electron densities (local/semilocal approximation)^^. Although the LDA and GGA 
functionals have proven to be rather universally applicable in theoretical materials science 
and achieve fairly good accuracy for ionization energies of atoms, dissociation energies of 
molecules, and cohesive energies, as well as bond lengths and geometries of molecules and 
solids^"—, this unexpectedly good performance for the ground-state properties of many 
materials is believed to be due to the partial error cancelation in the exchange and correlation 
energies22. In order to improve the performance of DFT on a theoretical basis, one approach 
is to add a portion of the non-local Hartree-Fock (HF) exchange to the local/semilocal 
approximation for exchange-correlation functional. In computational solid state physics, a 
breakthrough was achieved by Heyd, Scuseria, and Ernzerhof, who proposed the HSE03 
functional defined by22i~— : 

i?™ = \EriSR) + lEr'^iSR) + Er^(LR) + E^^. (1) 

The HSE03 functional benefits from the range separation of the HF exchange into a short- 
(SR) and long-range (LR) contribution and replaces the latter by the corresponding DFT 
exchange part. The IISE03 functional as well as its corrected form the HSE06 functional^^, 
have been extensively applied to calculate ground-state properties of solid and molecular 
systems, as well as adsorption energies for molecules, and have proven to be better than the 
LDA and GGA functionals in describing finite molecular systems. 

The results presented in this work are obtained using the projector augmented wave 
(PAW) method^ as implemented in the Vienna ab initio simulation package (VASP)^. For 
standard DFT calculations the xc energy is treated within the GGA using the parameter- 
ization of Perdew, Burke, and Ernzerhof (PBE)^i^ to compare with previously published 
calculations. In the present HSE calculations, the HSE06 hybrid functional^ is applied, 
where the range separation parameter is set to 0.2 A^^ for both the semilocal as well as the 
nonlocal part of the exchange functional. From now on, this particular functional will be 
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referred to as HSE. A plane wave basis set with a cutoff energy of 400 eV is adopted, and 
the scandium Ss^Sp^As^Sd^ and oxygen 2s^2p^ electrons are treated as fully relaxed valence 
electrons. A Fermi broadening^ of 0.05 eV is chosen to smear the occupation of the bands 
around the Fermi energy (Ej) by a finite-T Fermi function and extrapolating to T = K. 
The supercell containing the scandium oxide clusters is chosen to be without symmetry, and 
the cell size along each direction is larger than 15 A. A quasi- Newton algorithm is used to re- 
lax the scandium and oxygen ions for all scandium oxide clusters, with the force convergence 
criteria of 0.01 eV/A in PBE calculations and 0.03 eV/A in HSE calculations. 

III. RESULTS AND DISCUSSIONS 
A. Geometrical structures 

Firstly, the ScO and two isomers of Sc02 are shown in Fig.l (a) and (b) respectively. 
The optimized bond length of ScO is 1.677 A in PBE calculation, and 1.659 A in HSE 
calculation. The PBE result is in good agreement with former GGA calculations^^. The 
HSE06 result is near the ones obtained from the SCF^^ and hybrid BILYP calculations^!^. 
Both GGA and HSE06 calculations agree well with the experimental bond length (1.668A)^. 
Moreover, The vibrational frequency of ScO is 1056 cm~^ in HSE calculation, which is near 
the experimental value of 965 cm~^. 

For Sc02, the obtained ground state is the obtuse triangle structure, with the apex angle 
of 117.16°. It is interesting to note that these two Sc-0 bonds are not equivalent and the 
bond lengths are 1.687 and 1.928 Arespectively. In PBE calculations, the triangle is isosceles, 
with the apex angle of 126.19° and the ScO bond length of 1.781 A , in agreement with that 
have been found by using the DFT/PBEi^l and DFT/BPW91 methods^^. It is clear that 
the HSE06 calculation reduces the symmetry of the ground state. Both GGA and HSE06 
demonstrates the doublet is the most stable electronic state for the ground state of the Sc02 
cluster, in good agreement with the experimental result. The next stable state of Sc02 is 
the 0-Sc-O linear structure with 0.2 eV higher than the ground state. 

The low-lying isomers of Sc20m (m=l~4) are presented in Fig. 2 (a)-(d). The ground- 
state structure of SC2O cluster is also an obtuse triangle. The apex angle is 109.51°, and 
the two Sc-0 bonds are both 1.797 A long, indicating that the cluster has a point 
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(0.00 eV, M=2) 

(b) 




(0.00 eV, M=2) (0.20 eV, M=2) 

FIG. 1: (Color online). The low-energy structures of (a) ScO and (b) Sc02 clusters. Grey and red 
balls represent the scandium and oxygen atoms, respectively. The bond lengths are in units of A. 
The number in parenthesis is the relative energy (in eV) with respect to the corresponding ground 
state. Note that the spin multiplicity is also given in the parenthesis. 

symmetry. The ground state of SC2O is a triplet, which is in agreement with previous 
DFT studiesi^. The next stable structure of SC2O is linear, and has a triplet electronic 
state. The SC2O2 cluster adopts a distorted square as its ground-state structure. The Sc-0 
bond lengths are 1.892, 1.698, 1.892 and 1.698 A respectively. Since the PBE calculation 
prefers the rhombus with symmetry, our result is in accordance with the rule that the 
HSE06 functional tends to reduce the point group symmetry of some clusters obtained from 
the PBE calculations. The ground-state of SC2O2 is singlet, and the next stable state is 
the triplet state of the same structure. The other two structures of SC2O2 in Fig. 2(b) 
actually have much higher electronic free energies. For SC2O3, the ground-state structure 
is the singlet trigonal bi-pyramid structure shown in Fig. 2(c). This structure has a D^^h 
symmetry in PBE calculation, which disappears in HSE calculation. Until now, our PBE 
calculational results are in good agreements with the results reported by Wang et al}^ . The 
main modification of the HSE method, as we can see from the above description, is revealing 
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a lower symmetry for the ground-state structure. In the previous work by Wang et al, the 
ground-state structure of SC2O4 is the one with a molecular O2 adsorbing at one corner of 
the rhombus SC2O2 cluster, which is however, proven to be the next stable structure for 
SC2O4. The more stable structure we revealed for SC2O4 is the adsorbing structure with an 
O2 molecule adsorbing on top of the Sc-Sc bond, as shown in Fig. 2(d). The ground and 
next stable states of SC2O4 are both singlet. 

The SC3O cluster has been reported to have a singlet ground state, with the structure of 
an oxygen atom adsorbing on top of an equilateral scandium triangle^. Our result of SC3O 
is the same with previous, in both PBE and HSE calculations. The next stable structure 
of SC3O is the one with an oxygen atom adsorbing at the bottom of an isosceles scandium 
triangle, having a singlet electronic state and a free electronic energy of 0.16 eV higher than 
that of the ground-state of SC3O. Except for SC3O, other Sc30m clusters have seldom been 
studied before. Only a systematic study on the SC3O6 cluster is reported very recently^^. 
The ground-state of SC3O2 is the doublet state of the adsorbing structure of a scandium 
atom on the SC2O2 cluster. We also notice that there are several structure for the SC3O2 
cluster having close free electronic energies, which are thus shown in Fig. 3(b). They are the 
pentagon, bi-pyramid and bi-triangle structures respectively. The SC3O3 cluster adopts a 
hexagon structure with a fourfold spin multiplicity as its ground state. The smallest 0-Sc-O 
angle is 101.59° in PBE calculation, and 102.78° in HSE calculation. We also listed several 
other structures of SC3O3 in Fig. 3(c), especially the next stable doublet state of the folded 
rectangular structure, which is only 0.04 eV higher in free electronic energy. For SC3O4, the 
ground-state structure can be seen as the adsorption of an oxygen atom at the center of the 
hexagonal SC3O3 cluster, or the adsorption of three oxygen atoms at the three bottom sides 
of the pyramid SC3O cluster. The next stable structure of SC3O4, which is 0.42 eV higher 
in free electronic energy, can be seen as replacing one oxygen atom in the hexagonal SC3O3 
cluster with two oxygen atoms. 

As shown in Fig. 4(a), the ground-state structure for SC3O5 can be seen as adding an 
oxygen atom on top of the scandium atom in the ground-state structure of SC3O4. The Sc-0 
bond lengths become a little smaller than that in SC3O4. The next stable structure of SC3O5 
is the adsorption of 5 oxygen atoms on a SC3 linear chain, whose free electronic energy is 
however, much larger than the ground-state. The 3 lowest-energy structures of SC3O6 are 
the adsorption structures of 6 oxygen atoms on the SC3 triangle. In the ground state, the 
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oxygen atoms are separated from each other, while in the next and third stable structure, 
two of the six oxygen atoms bond together. We also find that the structures contain a Scs 
chain always have much larger free electronic energies. 

B. Important features from the structural studies 

Based on the above systematic study of Sc„Om clusters, we can see some features for 
both the oxidation pattern of scandium, and the influences of HSE calculation. 

For the oxidation of scandium, we conclude two important features. The first one is that 
the Sc„Om. cluster can be seen as adding m oxygen atoms to a Sc„ cluster, with the Sc^ 
cluster to be scandium dimer for n=2, and scandium triangle for n=3. The other feature 
is that the oxygen atoms in the ground-states of Sc„Om clusters are all bondless with each 
other, only one 0-0 bond forms in SC2O4. The fact that oxygen atoms do not bond with 
each other indicates that scandium oxide clusters are different from lead oxide clusters, in 
which oxygen atoms can form O2 or O3 units^. The structural features of Sc„Om clusters 
tell us that they can be built by adding m oxygen atoms separately to a Sc„ cluster. 

As for the improvements of the HSE method, the most important one is that the symmetry 
of the cluster is often lower in HSE calculation than in PBE calculation. The lengths of 
Sc-0 bonds always differ from each other in a cluster after geometry optimization using the 
HSE method. The second feature concluded from the above geometrical studies is that the 
HSE method does not change the relative stability between different structures. It means 
that geometry optimization using the PBE type functional will not yield wrong ground-state 
structures for scandium oxide clusters. At last, although the electronic structures in PBE 
and HSE calculations are the same for most Sc^O^ clusters, we do see differences in the 
ground-state electronic structures of SC3O2 and SC3O3. In PBE calculations, the magnetic 
moment is 3 for SC3O2, and 1 yU^ for SC3O3. And corresponding the magnetic moments 
obtained from HSE calculations are 1 //^ for SC3O2, and 3 /i^ for SC3O3. 

C. Fragmentation channels and dissociation energies 

The stability of scandium oxide clusters with different sizes and stoichiometrics is required 
to illustrate the growth pattern of various nanostructures and to understand even the oxida- 
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TABLE I: The most favorable fragmentation channels and dissociation energies {AE, in units of 
eV) of the Sc„Om (l<?^<3, 1 < m < 2n) clusters. All results are obtained by employing the 
HSE method. 



cluster 


Fragmentation channels 


AE 


ScO 


Sc+0 


+4.29 


ScOs 


ScO+0 


+1.41 


ScaO 


Sc+ScO 


+2.57 


SC2O2 


ScO+ScO 


+4.45 


SC2O3 


SC2O2+O 


+3.79 


SC2O4 


SC2O3+O 


+0.12 


SC3O 


SC+SC2O 


+1.78 


SC3O2 


SC+SC2O2 


+1.80 


SC3O3 


SCO+SC2O2 


+3.60 


SC3O4 


SCO+SC2O3 


+5.60 


SC3O5 


SC3O4+O 


+2.03 


ScsOe 


SC3O5+O 


+1.88 



tion behavior of the pristine scandium clusters. On the other hand, the study of the stability 
is helpful for finding the candidates of the building block of the cluster-assembled materials. 
Therefore, in the following, we evaluate the fragmentation energy of Sc„Om clusters. 

When a cluster A is dissociated into B and C fragments (i.e., A ^ B + C), the fragmen- 
tation energy is defined as AE = Eb + Ec — Ea, where the Eb, Eq, and Ea_ are the free 
electronic energies of clusters B, C, and A, respectively. In addition, the half of the total en- 
ergy of an oxygen molecule is adopted as the reference energy for atomic oxygen. Note that 
a fragmentation process is exothermic (endothermic) if the associated fragmentation energy 
is negative (positive). The dissociation energy AE, defined as the fragmentation energy of 
the most favorable channel, is calculated and presented in Table I. In general, the cluster 
with large positive AE has great stability, and that with small positive or even negative 
AE is not stable and tends to dissociate. Moreover, the frequently observed fragmentation 
products are believed to be stable^"^. Figure 5 shows the AE as a function of the oxygen 
atom numbers for for all clusters, revealing the underlying relationship between stability 
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and stoichiometry of these clusters. The essential features can be discussed as follows. 

Firstly, for the dissociation of Sc-rich clusters, the Sc atom is always one of the frag- 
ments. To make a further validation, we also study possible fragmentation ways for the 
SC4O cluster, using the HSE calculational method. It is found that the most energetically 
favored fragmentation channel is Sc40^Sc-|-Sc30, in which single Sc atom is also a frag- 
ment. Secondly, we do not see oxygen molecules in the fragments of Sc„0,„ clusters, the 
0-rich clusters favor the fragmentation channels that have a single oxygen atom as a prod- 
uct. Lastly, the dissociation energies of ScO, SC2O2, SC2O3, SC3O3, and SC3O4 clusters are 
larger than 3.60 eV. This indicates that the monoxide-like and sesquioxide-like scandium 
oxide clusters, and oxide clusters between them are remarkably stable. As a result, these 
clusters arc frequently observed in the fragmentation channels shown in Table 1. We also 
notice that the SC3O4 cluster has an enormous dissociation energy, and is more stable than 
the monoxide-like SC3O3 cluster. And at another side, the monoxide-like SC2O2 cluster, with 
a larger dissociation energy, is more stable than the the sesquioxide-like SC2O3 cluster. 

D. Electronic energy levels 

The electronic structures of transition metal oxide clusters are important characters for 
their experimental detections and chemical applications. Especially, their magnetic proper- 
ties are very important to reflect the behaviors of the (i-shell electrons. We thus calculate 
the ground-state electronic energy levels of Sc„0^ clusters, using the HSE method. Figure 
6 shows the obtained energy levels for different spin polarized electrons, in which spin-up 
and spin-down represent for the majority and minority spin states respectively. We find that 
the Sc„Om clusters with odd-number Sc atoms are all magnetic. For the Sc20m clusters, 
while the SC2O is magnetic, the other SC2O2, SC2O3, and SC2O4 clusters are nonmagnetic. 
The energy gap between the highest occupied molecular orbital (HOMO) and the lowest 
unoccupied molecular orbital (LUMO) is 1.82, 3.36, and 3.20 eV for the SC2O2, SC2O3, and 
SC2O4 cluster respectively. The large band gaps indicate that the SC2O2, SC2O3 and SC2O4 
clusters are chemically very stable. In contrast, although the SC3O4 cluster is energetically 
very stable, its energy gap of spin-up electrons is as small as 0.49 eV, indicating that it is 
chemically active, and ready to interact with other particles or materials. 

To systematically investigate the electronic structures of Sc„0^ clusters, we then calcu- 
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late the projected density of states (PDOS) for the ScnO^ clusters, which are spin-resolved. 
Figures 7(a)-7(l) hst the obtained s-, p-, and d-PDOS of oxygen and scandium atoms, for 
the Sc„0^ clusters respectively. We can clearly see that for all the studied Sc„Om clusters, 
there are strong hybridizations between oxygen 2p and scandium 3d electronic states, and 
they contribute most of the electronic states around the Fermi energies. For Sc-rich clusters 
like SC2O, SC3O, and SC3O2, the HOMO and LUMO are both composed of Sc-3d electrons. 
For the three nonmagnetic SC2O2, SC2O3, and SC2O4 clusters, the HOMO and LUMO are 
contributed by 0-2p electrons and Sc-3d electrons, respectively. And for 0-rich clusters hke 
SCO2, SC3O5, and SC3O6, the HOMO and LUMO are both 2p electronic states of oxygen. The 
electronic structure of ScO is, however, different from the other clusters, for its HOMO of 
spin-up electrons and LUMO of spin-down electrons are contributed by 4s electronic states 
of scandium. 

Since the obtained ground-state spin configurations are different for the SC3O2 and SC3O3 
clusters in PBE and HSE calculations, we draw their PDOS together to compare the differ- 
ences in electronic-state descriptions of PBE and HSE methods, which are shown in Figs. 
8(a)-8(d). In the p-d hybridization area from -7.0 eV to -3.0 eV below the Fermi energies, we 
can sec that the hybridized peaks arc more separate in HSE calculations for both the SC3O2 
and SC3O3 clusters. The electronic states around the Fermi energy are both Sc-3d states for 
the two chosen clusters. One can see from Fig. 8 that there are always more localized Sc-3d 
peaks around the Fermi energy in HSE calculations than in PBE calculations. It means that 
the HSE method can lead to more locahzed descriptions for Sc-3d electronic states in Sc^O^ 
clusters. Considering that the PBE and other standard GGA methods rely on the xc energy 
of the uniform electron gas, and thus are expected to be useful only for systems with slowly 
varying electron densities, we think that the HSE descriptions on the Sc-3d states are more 
reasonable. 

IV. CONCLUSIONS 

In this work, first-principles calculations with the HSE method have been performed to 

study the geometries, stabilities and electronic structures of small Sc„Om clusters (n=l- 
3,m=l-2n). Based on an extensive search, it is found that the lowest-energy structures 
of all these clusters can be obtained by the sequential oxidation of small "core" scandium 
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clusters, with the adsorbing oxygen atoms separating from each other. 

The fragmentation analysis reveals that the ScO, SC2O2, SC2O3, SC3O3, and SC3O4 clusters 
have great stability. This suggests that these clusters might be used as candidates of the 
building block of cluster-assembled materials. The fragmentation of Sc-rich clusters is found 
to include a scandium atom as a product, while the fragmentation of 0-rich clusters is found 
to include an oxygen atom as a product. Besides, the above four extremely stable clusters 
can also be frequently seen in the fragmentation products of other ScnO^ clusters. 

Through electronic structure calculations and wavef unction analysis, we reveal that most 
Sc oxide clusters have magnetic ground state except SC2O2, SC2O3, and SC2O4 clusters. 
The HOMO and LUMO of the three nonmagnetic clusters arc composed of 0-2p and Sc- 
3d electronic states, respectively. For the Sc-rich clusters, the HOMO and LUMO are 
contributed both by Sc-3d electrons, while for 0-rich clusters, the HOMO and LUMO are 
contributed by 0-2p electrons. 

In comparison with standard PBE calculations, the HSE method is superior because 
it can correct the wrong symmetries and electronic configurations in PBE results of some 
clusters. 
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FIG. 2: (Color online). The low-energy structures of (a) SC2O, (b) SC2O2, (c) SC2O3, and (d) SC2O4 
clusters. Grey and red balls represent the scandium and oxygen atoms, respectively. The bond 
lengths are in units of A. The number in parenthesis is the relative energy (in eV) with respect to 
the corresponding ground state. Note that the spin multiplicity is also given in the parenthesis. 
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FIG. 3: (Color online). The low-energy structures of (a) SC3O, (b) SC3O2, (c) SC3O3, and (d) SC3O4 
clusters. Grey and red balls represent the scandium and oxygen atoms, respectively. The bond 
lengths are in units of A. The number in parenthesis is the relative energy (in eV) with respect to 
the corresponding ground state. Note that the spin multiplicity is also given in the parenthesis. 
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FIG. 4: (Color online). The low-energy structures of (a) SC3O5 and (b) ScaOg clusters. Grey and 
red balls represent the scandium and oxygen atoms, respectively. The bond lengths are in units 
of A. The number in parenthesis is the relative energy (in eV) with respect to the corresponding 
ground state. Note that the spin multiplicity is also given in the parenthesis. 
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FIG. 5: (Color online). Dissociation energies of the Sc„Om clusters as a function of the oxygen 
atom number m. The corresponding fragmentation channels are listed in Table I. 
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FIG. 6: (Color online). The electronic energy levels for (a) ScO, (b) Sc02, (c) SC2O, (d) SC2O2, 
(e) SC2O3, (f) SC2O4, (g) SC3O, (h) SC3O2, (i) SC3O3, (j) SC3O4, (k) SC3O5, and (1) ScsOe clusters, 
obtained from HSE calculations. Spin-up and spin-down electronic states are denoted by black and 
red lines, respectively. The Fermi levels are denoted by the blue dotted lines. 
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FIG. 7: (Color online). The spin-resolved projected density of states for the (a) ScO, (b) Sc02, (c) 
SC2O, (d) SC2O2, (e) SC2O3, (f) SC2O4, (g) SC3O, (h) SC3O2, (i) SC3O3, (j) SC3O4, (k) SC3O5, and 
(1) SC3O6 clusters. The Fermi energies are all set to be zero. The electronic states of O and Sc are 
shown in dotted and solid line respectively, while s, p and d electronic states are shown as black, 
blue and red lines respectively. 
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FIG. 8: (Color online). The spin-resolved projected density of states for the SC3O2 cluster in (a) 
PBE and (b) HSE calculations, and the SC3O3 cluster in (c) PBE and (d) HSE calculations. The 
Fermi energies are all set to be zero. The electronic states of O and Sc are shown in dotted and 
solid line respectively, while s, p and d electronic states are shown as black, blue and red lines 
respectively. 
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